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, We analyze the superfield constraints of the Z)=4, iV=3 SYM-theory using light-cone gauge 

conditions. The SU(3)/U(1) x U(l) harmonic variables are interpreted as auxiliary spectral 
parameters, and the transform to the harmonic-superspace representation is considered. Our 
nilpotent gauge for the basic harmonic superfield simplifies the SYM-equations of motion. A 
partial Grassmann decomposition of these equations yields the solvable linear system of iterative 
equations. 
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% ; 1 Introduction 

It is well known that the geometric superspace constraints of the D = 4, N = 3 SYM- 
theory are equivalent to the equations of motion for the component fields [1] . The special 
projections of these superfield constraints can be interpreted as conditions of zero curva- 
tures assosiated with the Grassmann covariant derivatives [2,3]. The possible connection 
of these zero-curvature conditions with integrability or solvability of the super- Yang-Mills 
(SYM) theories with 16 or 12 supercharges has been discussed more than 20 years in 
the framework of different superfield approaches (see, e.g. [2-6]). Here we shall use the 
harmonic-superspace approach [3] to analyze the solutions of the N = 3 SYM-equations. 

The harmonic-superspace (HSS) method has been introduced first for the solution of 
the D = 4, N = 2 off-shell superfield constraints [7]. Harmonic variables are analogous 
to twistor or auxiliary spectral variables used in integrable models. Harmonic and twistor 
methods give the explicit constructions of the general solutions to zero-curvature condi- 
tions in terms of independent functions on special analytic (super) spaces satisfying the 
generalized Cauchy-Riemann analyticity conditions. For instance, the off-shell N = 2 
superfields satisfy the conditions of the Grassmann (G-) analyticity. 

In the standard harmonic formulation of the D = 4, N = 2 SYM-theory, the basic 
harmonic connection is G-analytic [7], and the 2-nd one ( via the zero-curvature condition) 
appears to be a nonlinear function of the basic connection. The N = 2 equation of motion 
is linearly dependent on the 2-nd harmonic connection, but it is the nonlinear equation 
for the basic connection [8]. It has been shown in Ref. [5] that one can alternatively choose 
the 2-nd harmonic connection as a basic superfield, so that the dynamical G-analyticity 
condition (or the Grassmann-harmonic zero-curvature condition) for the first connection 
becomes a new equation of motion. We shall use below a similar change of the basic HSS 
variables for the N = 3 SYM-theory. 

In the HSS-approach to the D = 4, N = 3 supersymmetry, the SU(3)/U(1) x £7(1) 
harmonics have been used for the covariant reduction of the spinor coordinates and deriva- 
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tives and for the off-shell description of the SYM-theory in terms of the corresponding 
G-analytic superfields [3]. The N = 3 SYM-equations in the ordinary superspace have 
been transformed to the zero-curvature equations for the harmonic gauge connections, 
however, nobody tried earlier to obtain the general solution of these equations. 

The alternative SL(2, C)-harmonic interpretation of the N = 3 SYM-equations and 
the corresponding harmonic zero-curvature equation for two G-analytic connections has 
been considered in Ref.[6]. We shall not discuss in this paper the problem of solving the 
equations of the SL(2, 0)-harmonic approach. 

An analogy with the self-dual SYM-solution [9] has been used recently for the analysis 
of the 10.D SYM-solutions [10]. These authors have considered the very useful light-cone 
gauge conditions which break the space-time symmetry of the SYM-equations. We shall 
consider the analogous light-cone gauge conditions which preserve the SU (3)-invariance 
and allow us to use the harmonic approach. 

It should be underlined that the N = 3 harmonic equations of motion contain three 
analytic connections, so they are more complicated than the SU (2)-harmonic equations 
using the single analytic connection. If one does not treat the 577(3)-harmonic connections 
as the basic variables, then the harmonic zero-curvature equations can be readily solved in 
terms of the non-analytic superfield matrix v ( bridge). In this approach, the G-analyticity 
conditions for the composed harmonic connections can be considered as the dynamical 
equations for this independent superfield v. 

It will be shown that the light-cone gauge conditions for the initial superfield con- 
nections correspond to the light-cone analyticity of the bridge matrix. We choose the 
nilpotent gauge for v and construct the nilpotent analytic harmonic connections in this 
representation. We demonstrate that the 1-st order harmonic bridge equations with the 
nilpotent connections produce also linear 2-nd order differential constraints for the basic 
matrices. 

This talk is based on our work [11]. 

2 Reduced-symmetry representation of D = 4, N = 3 
SYM constraints 

The covariant coordinates of the D — 4, N = 3 superspace are 

Z M = Qf, 6™) , (1) 

where a, a are the SL(2, C) indices and i — 1, 2, 3 are indices of the fundamental repre- 
sentations of the group SU(3). 

We shall study solutions of the SYM-equations using the non-covariant representation 
and the new notation for these coordinates 

,,3 = 22 ± 3 12 1,-2 - 2i 1 -2 

t + x , x = x — t — x , y = x = x + %x , y = x = x — tx , 
, {6 l+ , = Q i(x . (2) 

based on the reduced symmetry 5*0(1, 1) x SU(3). 

The 5*0(1, 1) weights (helicities) of these coordinates are 

w(x + ) = 2 , w(x=) = -2 , w(y) = w(y) = , w(6f) = w(^ ± ) = ±1 . (3) 
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The reduced representation of the algebra of D = 4, N = 3 spinor derivatives is 



{D k + , D l + } — , {D k+ , D l+ } = , {£>* , A+} = 2i«J* 9 4 , 

DL} = , {A-, A-} = , , A-} = 2i«y*9 = , (4) 

A-} = 2i6? d y , {D k _, A+} = 2z5fa, , D l _} = {D k+ , A-} = . 

The last two relations can contain, in principle, six central charges, however, we shall 
consider the basic superspace without central charges. The general N = 3 superspace has 
the odd dimension (6,6) in our notation with the reduced SO(l, 1) symmetry. 

Let us consider the (4|6, 6)-dimensional superspace gauge connections A(z) and the 
corresponding covariant derivatives V 

V l ± = D { ± + A\ , V l± = A± + A i± , (5) 
V 4 = d + + A + , V= = d= + A= , V y = ^ + A y , = d y + A y . 

The D = 4, iV = 3 SYM-constraints [1] have the following reduced-superspace form: 

{V*,V^} = 0, {V fc+ ,V i+ } = 0, {V*,V {+ } = 2i#V + , (6) 

{V* , V'_} = W w , {V* , V,_} = 2itf V„ , (7) 

{V* , V i+ } = 2i#V„ , {V fe+ , V/_} = W w , (8) 

{V* , V!_} = , {V fc _, V,_} = , {V fc _, V,_} = 2z#V = , (9) 



where Wki and are the gauge-covariant superfields constructed from the gauge connec- 
tions. This reduced form of the 4D constraints is convenient for the study of dimensionally 
reduced SYM-subsystems in D < 4. 

Let us analyze first Eqs.(Q) combined with the relations 

[V*, V + ] = [V fc+ , V + ] = . (10) 

The simplest light-cone gauge condition is 

A k + = 0, A k+ = 0, A^=0. (11) 

The analogous gauge conditions excluding the part of connections has been considered in 
Ref.[9] for the self-dual AD SYM-theory and in Ref.flO] for the 10D SYM equations. 

Stress that the light-cone gauge preserves the S'?7(3)-symmetry of the N = 3 superfield 
constraints, so we can use the harmonic-superspace method for non-trivial equations 
in this gauge. 



3 Harmonic-superspace equations for the nilpotent 
bridge matrix 

The SU(3)/U(1) x £7(1) harmonic superspace has been introduced in Ref. [3] for the off- 
shell description of the AD N = 3 SYM-theory. The dynamical equations in this approach 
have been transformed into the set of pure harmonic equations for the G-analytic super- 
field prepotentials. It should be stressed that nobody tried earlier to analyze in details 
the solutions of these harmonic equations. 
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Now we shall study the reduced-symmetry version of the SU(3)/U(1) x U(l) harmonic 
superspace which allows us to consider the non-covariant gauges and the dimensional 
reduction. 

The SU(3)/U(1) x U(l) harmonics [3] parameterize the corresponding coset space. 
They form an SU(3) matrix u\ and are defined modulo U(l) x U(l) transformations 

«J=t4 1,0 \ u^u^, u 3 = u!- '- 1 ) , (12) 

where the lower index i describes the triplet representation of SU (3), and the upper indices 
1, 2 and 3 correspond to different combinations of charges. The complex conjugated 
harmonics have opposite U(l) charges and reverse positions of indices 

u\ = v,*- 1 ® , 4 = u^ 1 '" 1 ) , ul = v*<M . (13) 

These harmonics satisfy the following relations: 



,jfcL,l„.2„,3 



e M u\uiui = 1 . (14) 
The SU (3)-invariant harmonic derivatives act on the harmonics 

d'juf = 8«ul , d^K = SWj , 

[9 I j,d£] = 8?dL-5 I L dK. (15) 
We shall use the special S'[/(3)-covariant conjugation 

uj = u\ , ul = u[ , u 2 = -u\ . (16) 

We can define the real analytic harmonic superspace H(4, 6|4, 4) with 6 coset harmonic 
dimensions and the following left and right coordinates: 

C = x= , y , y , ef , ef , e l± , e 2± ) , 

X*= x++ i(9+9 3+ - 9+9 1+ ) , X= = x= + i(9 3 9 3 ~ - flf/0 1 -) , 

y = y + t(e+e 3 - - e+e 1 -) , y = y + i(e^e 3+ - e^e 1+ ) , (17) 

where 9f = 9^u) , 9 ±J = 9 ±k u{. 

The CR-structure in if(4, 6|4, 4) involves the G-derivatives 

Dl, D 3± , (18) 

which commute with the harmonic derivatives D\, Dl and D\. 

These derivatives have the following explicit form in the analytic coordinates: 

£>£ 1,0) = Dl = dl = d/dOf D { ±' 1] = D 3± = 8 3 ± = d/d9 3± , (19) 

z^ 2 '- 1 ) = d\ = dl + ie+P+d^+ie+p-dy + ie 2 e 1+ d Y + ie 2 e 1 -d = 
-e+d l + - e 2 d l _ + e l+ B 2+ + e 1 -d 2 - , 

D*" 1 ' 2 ) = Dl = dl + i9+9 2+ d + + i9+9 2 ~d Y + i9 3 9 2+ d Y + i9 3 9 2 -& = 
-9+d\ - 9 3 d 2 _ + 9 2+ d 3+ + 9 2 ~d 3 . , (20) 
D^ l) = D\ = dl + 2i9+9 1+ 8 + + 2i9+9 1 ~d Y + 2i9 3 9 1+ B Y + 2i9 3 9 1 ~d = 
-9+dl - 9 3 d l _ + 9 l+ B 3+ + p-Ba- , 



where d^=d/dX* d = = d/dX=, d Y = d/dY and d Y = d/dY. 

It is crucial that we start from the light-cone gauge conditions ( |TT| ) for the N = 3 
SYM- connect ions which break SL(2, C) but preserve the S'f/(3)-invariance. Consider 
the harmonic transform of the covariant Grassmann derivatives in this gauge using the 
projections on the S[/(3)-harmonics 

V 7 + =u\D\ = D\, V/+ = u)D i+ = D I+ , {D J + , D K+ } = 2i5 I K d + , (21) 
V T _ = ulVi = D T _ + AL, Vj- = u\Vi- = Di- + Ai- , (22) 

where the harmonized Grassmann connections A[_ and Ai- are defined. 

The S't/(3)-harmonic projections of the superfield constraints (^-§) can be derived from 
the basic set of the N = 3 G-integrability conditions for two components of the harmonized 
connection: 

D\A l _ = D 3+ A\ = D\A^ = D 3+ A 3 ^ = (23) 
D\A\ + (Alf = , D 3 _A 3 - + (^3-) 2 = , 

D\A 3 - + D 3 _A l _ + {Al, A 3 ~} = . (24) 

All projections of the SYM-equations can be obtained by the action of harmonic deriva- 
tives D T K on these basic conditions. 

These Grassmann zero-curvature equations have the very simple general solution 

A\[v) = e- v D\e v , A 3 -{v) = e~ v D 3 _e v , (25) 

where the bridge v is the matrix in the Lie algebra of the gauge group. This superfield 
matrix satisfies the additional constraint 

(Di, D 3+ )v = , (26) 



which is compatible with the light-cone representation fl21|). Thus, v does not depend on 
the Grassmann coordinates Of and # 3+ . 

Consider the gauge transformations of the bridge 

e v => e x e v e Tr , (27) 

where A G H(4, 6|4, 4) is the (4,4)-analytic matrix parameter, and the parameter r r does 
not depend on harmonics. The matrix e v realizes the harmonic transform of the gauge 
superfields A±, A k ± in the central basis (CB) to the equivalent set of harmonic gauge 
superfields in the analytic basis (AB). We have partially fixed the CB-gauge invariance 
in fl2"E|), and the A-gauge transformations of v will be used below. 

The dynamical SYM-equations in the bridge representation fl25|) are reduced to the 
following harmonic differential conditions for the basic Grassmann connections: 

il r>2 nl\ ( /il 



(Dl Dl Dl)(A\(v), A 3 4v))=0. (28) 



Using the off-shell (4, 4)-analytic A-transformations one can choose the non-supersym- 
metric nilpotent gauge condition for the superfield bridge 

v = 0-b 1 + p-b 3 + e^P'dl , v 2 = 07/0 3 ~ [b 3 , b 1 ], v 3 = , (29) 

e~ v — i — v + ^v 2 — i — e^b 1 - p-b 3 + e^-{hp 3 , b 1 ] - d\) , (30) 
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where the fermionic matrices b l ,b 3 and the bosonic matrix d\ are analytic functions of 
the coordinates (. 

Note that the nilpotent gauges for the harmonic bridges and connections are possible 
for the harmonic formalisms with the off- shell analytic gauge groups only [3,7]. 
The conditions for the S77(n)-bridge are 

Tr v = , v* — —v . (31) 

The matrices fe 1 ,^ and d\ have the following properties in the gauge group SU{n): 

Tr b 1 = , Tr b 3 = , Tr d\ = , (32) 
(&*)* = 6s, (4) f = "4- (33) 

Let us parameterize the Grassmann connection A l _{v) and A 3 - in terms of the basic 
analytic matrices ft 1 , 63 and d\ ( |29|) 



^L(«) = e-Die B = b 1 - e^b 1 ) 2 + 9*-fl + e^'[b\ fl] , (34) 
A 3 - = e- v D 3 _e v = b 3 + d^fl - 6 3 -(b 3 ) 2 + B^~\JIM , (35) 

where the following auxiliary superfields are introduced: 

fl = d\-\{b\h}, fl = -d\- l -{b\h} . (36) 

Equations (D\, D 2 )A l _{v) = generates the following independent relations for the 
(4,4)-analytic matrices: 

D]b l = -d^b 1 ) 2 , (37) 
Dlb 1 = -9 2 -fl , (38) 
D\fl = O^flb 1 ] , D 2 3 fl = 0. (39) 

Equations (D\,Df)A 3 -{v) = are equivalent to the relations 

Dfb 3 = 6 2 f^ I)% O 1 [b 3 f . (40) 

£2/3=0, D 2 Jl = e 2 -\l 3 Jl}. (4i) 

In the case of gauge group SU(n), the last equations are not independent, they can be 
obtained by conjugation from the equations for b 1 and /g. 

It is useful to derive the following relations for the matrices b 1 and b 3 which do not 
contain the auxiliary matrices d\, f 3 or /3 1 : 

e 2 D\{b\ 63) = , e 2 -Dl{b\ 63) = , (42) 

e^Dlb 1 + 6 2 -Dlb 3 = 9 2 6 2 -{b\ b 3 } . (43) 
Solutions of the linear equations for matrices /3 1 and satisfy the subsidiary condition 

+ ft = -{b\h} ■ (44) 

It is important that all these equations contain the nilpotent elements 6 2 or 9 2 ~ in the 
nonlinear parts, so they can be reduced to the set of linear iterative equations using the 
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partial Grassmann decomposition. In particular, the nilpotency of the basic equations 
yelds the subsidiary linear conditions for the coefficient functions 

D 1 2 D 1 2 (b\b 3 ,dl)=0, D 2 3 Dl(b\b 3 ,dl) = 0. (45) 

The harmonic linear equations for the analytic superfields fa 1 , 63, d\ have simple (short) 
solutions with the finite number of the harmonic on-shell field components f\. This har- 
monic shortness is an important restriction on the structure of the SYM-solutions. 



4 Analytic representation of solutions 

Remember that the following covariant Grassmann derivatives are flat in the analytic 
representation of the gauge group before the gauge fixing: 

e v Vle- v =Vi = Di, e v V 3± e~ v = V 3 ± = D 3 ± (46) 

The harmonic transform of the covariant derivatives using the matrix e v fl2"5| ) determines 
the composed on-shell harmonic AB-connections 

V J K = e v D' K e-" = D I K + V£(v) , 

Vk(v)=e v (D I K e~ v ) . (47) 

Note that the harmonic connections in the bridge representations satisfy automatically 
the harmonic zero- curvature equations. 

It is evident that the basic equations (|28[) are equivalent to the following set of the 
dynamic G-analyticity relations for the composed harmonic connections: 

pi, A.-) (vj(y), V 3 2 (v), Vj(vj) = , (48) 

The positive-helicity analyticity conditions are satisfied automatically for the bridge in 
the gauge (p6|) . 

The dynamical G-analyticity equation fllcf) in gauge ( p9f) is equivalent to the following 
harmonic differential bridge equation: 

V 2 \v) = 2 b l , (Vi) 2 = (49) 
Die-* = e~ v Vi . (50) 

where the manifestly analytic nilpotent representation for the harmonic connection is 
used. 

The 2-nd on-shell harmonic connection is also nilpotent V 3 (v ) = —6 2 ~b 3 . 
The spinor AB-connections can be calculated via the non-analytic harmonic connections 
by analogy with the N = 2 formalism of Ref . [8] 

a% = -DlV 2 , al = -D 1 ± V 1 3 , (51) 
a 2± = D 3 ±V 2 3 , oi± = D 3± V? (52) 

where V 2 , V 3 and V 2 are the non-analytic harmonic connections. 

2 The analogous short harmonic N — 3,4 Abelian superfields have been considered in Ref. [13]. 
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The 3-rd harmonic analytic connection can be readily calculated 

Vj = D\Vl - D 2 3 Vi + [V 2 \ V 2 ] = e.b 1 - B^h , (53) 

where Eq.(^) is used. 

The harmonic connection V 2 can be written in terms of the superfield b l only 

V 2 = -O^Dfb 1 . (54) 

The conjugated connection V 2 = — {V^Y contains matrix 63 only. 
These connections satisfy the partial G-analyticity conditions 

D 3± V 2 = , DiV 2 3 = . (55) 

It is convenient to define the AB-superfield strengthes 

= -D\D\V 2 = -D 2 + b l , (56) 

where Eq.(Q) is used. 

Stress that the single coefficient matrix b 1 generates the family of the AB-geometric 
objects: V^, V 2 , a\ and w 12 . The conjugated 6 3 -family of superfields contains V^, V 2 3 , a 2 ± 
and W23. 

The superfield w 12 satisfy the (4,2)-dimensional G-analyticity conditions 

D l ± w 12 = D 3± w 12 = D 2 ± w 12 + [4, w 12 ] = (57) 
and the non-Abelian H-analyticity conditions 

(V 1 2 ,VlV 1 3 )w 12 = Q. (58) 

The (4,2)-analytic superspaces have been considered earlier in Refs. [12,13]. 

It should be noted that function d\ (or /3 1 ) is an auxiliary quantity in the framework 
of the analytic basis, since all harmonic and spinor connections and tensors of this basis 
can be expressed in terms of superfields b 1 and 63 only. Nevertheless, the construction of 
d\ is important for the transition to the central basis. 

Let us analyze the analytic equations for the basic fermionic (4,4) matrices b 1 and 63. 
The nonlinear terms in these equations contain the negative-helicity Grassmann coordi- 
nates 6 2 ,6 3 ,9 l ~ and 6 2 ~ , so the partial decomposition in terms of these coordinates is 
very useful for the iterative analysis of solutions. Equations for the auxiliary matrix d\ 
do not give additional restrictions on b 1 and 63. 

Consider first the decomposition of the harmonic derivatives and define the harmonic 
derivatives on the (4,0) analytic functions depending on the analytic Grassmann variables 
6+,6t,6 1+ ,6 2+ 



bl 


- d 1 


+ i9+6 1+ d^-6+dl^ 


- o 1+ d 2+ , 






= dl 




- 2+ B 3+ , 






= d\ 


+ 2i6+6 1+ d^-6+d 1 + 


+ 1+ B 3+ . 


(59) 
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The (4,0) decomposition of the (4,4) analytic matrix functions has the following form: 



b l = p 1 + e 2 B 12 + 9 3 b 13 + p-b + e 2 ~B\ + 9 2 9 3 p° + 9 2 9 l ~p 2 + 9 3 9 2 ~p\ 3 

+9 3 ~9 1 -p 3 + 9 1 ~9 2 ~p 2 + 2 6 2 -ri l + Q 2 Q^¥-B 2Z + Q^¥-Q 2 -B\ + 9 2 9 3 9 2 -C 13 

+6 2 e 1 -0 2 -c° + e 2 e^e l -e 2 -rf . (60) 

The analogous decompositions can be written for 63 and d\. 

It is easily to show that a part of the (4,0) coefficients can be constructed as the algebraic 
functions of the basic set of independent (4,0) matrices, for instance, 

B\ = -b\B° - i9+d Y p 1 , (61) 

B 12 = D 2 3 B 13 + i6 2+ B Y f3 l , (62) 

p 2 = t) 2 3 p 3 - i9 2+ d Y B° , (63) 

r! 1 = -Dip 2 + i9 1+ B Y B° + i9+d Y B 12 + \J3\ B°] . (64) 

The independent (4,0) matrix functions are 

5°, fi 13 , B 23 , (C° - C°), B 3 , p\ fa P 3 , p\ V 3 . (65) 
The (4,0) matrix of the dimension I — —1/2 satisfies the linear equations 

DIP 1 = Djp 1 = . (66) 
The equations for the I — — 1 matrices B 13 and B° are 

D\B 13 = , {Dl) 2 B 13 = , 

D\B 13 = -2i9 1+ B Y p x - (p 1 ) 2 , (67) 
D 2 3 B° = , (bl) 2 B° = , 

bl(B°-B°) =t9+d Y p 1 + t9 1+ d Y p 3 + {P\p 3 } . (68) 

The inhomogeneous linear equations for B 13 and B° contain sources with the functions 
P 1 and Ps = (P 1 )^ calculated on the previous stage. 

The iterative equations for the (4,0) matrices with / < — 1 can be analyzed analogously. 
Each independent iterative equation is manifestly resolved in terms of the harmonic deriva- 
tives of the corresponding function and the sources of these equations can be calculated 
on the previous stage of iteration. 

Thus, it can be shown easily that the basic iV = 3 harmonic equations for the (4,4) 
analytic moduli functions with the nilpotent nonlinear terms are equivalent to the fi- 
nite number of the linear iterative (4,0) equations which contain non-Abelian sources 
constructed from the solutions of the previous step of iteration. 

Note that all iterative equations are simplified essentially for the two-dimensional solu- 
tions which do not depend on variables Y and Y. 



5 Conclusions and acknowledgment 

We have described the harmonic transform of the iV = 3 SYM-equations of motion in 
the standard superspace to the differential HSS equations for the bridge matrix e v which 
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connects different representations of the gauge superfields. The light-cone nilpotent gauge 
condition for matrix v simplifies significantly the analysis of the bridge equations. This 
condition yields the nilpotent harmonic-analytic gauge connections and the corresponding 
linear 2-nd order differential conditions for the basic matrices. The nilpotency of nonlinear 
terms in the basic HSS equations allow us to consider the simple iterative procedure based 
on the partial decomposition in Grassmann variables of the negative helicity. The finite 
set of these solvable linear iterative equations can be used for the explicit construction of 
the N = 3 SYM-solutions in the harmonic superspace and the on-shell gauge superfields 
in the ordinary superspace. 

Let us discuss shortly the problems of solving equations of the N — 0, 1 and 2 subsys- 
tems of the N = 3 theory. One can study the problem of constructing solutions with the 
reduced number of fermion or scalar fields using the explicit constructions of the N = 3 
solutions, although the formal reduction of the stable and regular N = 3 solution could, 
in principle, correspond to unstable or irregular solutions with lower supersymmetry. 

The N = 2 SYM-constraints together with equations of motion have the Grassmann- 
harmonic zero-curvature representation in the framework of the S , t/(2)-harmonic super- 
space [5]. It should be underline that the nilpotent bridge representation similar to (|29|) 
is also very useful for solving the N = 2 SYM-equations. The corresponding N = 2 
solutions will be discussed elsewhere. 

The author is grateful to J. Niederle for collaboration and to E. Ivanov and E. Sokatchev 
for discussions. 

This work is supported by the Votruba-Blokhintsev programme in Joint Institute for 
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and NATO-PST.CLG-974874. 
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